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1.  Introduction. 


Let  X,  Y  be  Banach  spaces  and  A  C  R"  a  bounded  interval.  Let  /’:AxX— ^ybea  smooth 
operator.  The  nonlinear  equation 

(1.1)  F(A,«)  =  0, 

with  parameters  A  €  A  is  called  parametrized  nonlinear  equations. 

Let  (A,  u)  €  A  X  A*  be  a  solution  of  (1.1).  Intuitively,  the  set  of  the  solutions  of  (1.1)  would 
form  n-dimensional  hypercurves  in  the  Banach  space  R"  x  X.  If  DuF{\,u)  €  C(X,Y),  the 
Frechet  derivative  of  F  with  respect  to  u,  is  an  isomorphism,  then,  by  the  implicit  function 
theorem,  the  above  intuition  is  correct,  i.e.  there  exists  a  locally  unique  branch  of  solutions 
around  (A, u),  and  the  branch  is  parametrized  by  A.  Such  branches  on  which  DuF(\,u)  is 
isomorphism  at  each  (A,u)  are  called  regular  branches. 

However,  if  Duf'(A,u)  is  not  an  isomorphism,  the  state  of  equilibrium  defined  by  (1.1) 
becomes  unstable  and  the  behavior  of  the  solutions  is  unpredictable;  the  hypercurve  of  the 
solutions  might  be  a  fold,  or  there  might  be  several  hypercurves  of  solutions  intersecting  at  that 
point.  The  folding  points  are  called  turning  points.  The  points  at  which  the  hypercurves  of 
solutions  are  intersecting  are  called  bifurcation  points.  (Note  that  the  definition  of  bifurcation 
points  given  by  some  authors  includes  turning  points.) 

In  this  paper  we  deal  with  the  parametrized  nonlinear  equation  F  :  A.  x  Hq{J)  — 
with  one  parameter  A  €  A  defined  by 

(1.2)  F’(A,«)  =  0,  (A,«)eAxifoV), 

(1.3)  <  F{X,u),v  >:=  +  f{X,x,u{x))v]dx,  Vu  6 

where  J  :=  (6,c)  C  R  is  a  bounded  interval,  and  a,  /  :  A  x  7  x  R  — *  R  are  sufficiently  smooth 
functions.  Since  F  is  a  second  order  differential  operator  in  divergence  form,  finite  element 
solutions  of  (1.2)  are  defined  in  a  natural  way. 

Brezzi,  Rappaz,  and  Raviart  [BBRl-3]  presented  a  comprehensive  work  on  the  numerical 
analysis  of  parametrized  nonlinear  equations.  They  first  proved  an  extended  implicit  function 
theorem  with  error  estimates  on  Banach  spaces.  Then,  using  the  implicit  function  theorem, 
they  obtained  several  results  of  a  priori  error  estimates  of  finite  element  solutions  [BRRl,2].  In 
[BRR3],  they  considered  approximation  of  solution  branches  around  bifurcation  points,  which 
will  not  be  dealt  with  in  this  paper. 
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Following  [BRR1-3|,  Fink  and  Rheinboldt  released  several  papers  about  numerical  analysis 
of  parametrized  nonlinear  equations  ([FR1,2],  [R],  and  references  therein).  While  the  formulation 
of  [BRRl-3]  was  rather  restrictive,  Fink  and  Rheinboldt  developed  their  theory  of  a  priori  error 
estimates  of  numerical  solutions  in  a  very  general  setting  using  the  theory  of  differential  geometry. 

Fink  and  Rheinboldt  employed  the  theory  of  Fredholm  operators.  Let  X  and  Y  be 
Banach  spaces  and  F  :  X  —  Y  3^  differentiable  mapping.  Then,  F  is  called  Fredholm  on  an 
open  sei  U  C  X  if  the  Frechet  derivative  DF(x)  satisfies  the  following  conditions  at  any  x  €  U: 

(1)  dimKerDF  is  finite, 

(2)  ImDF  is  closed, 

(3)  dimCokeri?F  is  finite. 

We  must  note  that,  in  the  above  prior  works  by  Brezzi,  et  al.  and  Fink-Rheinboldt,  only 
mildly  non/inear  equations  were  considered.  If  a(A,i,y)  in  (1.3)  is  nonlinear  with  respect  to 
y,  the  operator  F  is  called  strongly  nonlinear  (quasilinear),  otherwise  it  is  called  mildly 
nonlinear  (semilinear). 

Following  the  above  prior  works,  we  here  develop  a  thorough  theory  of  a  priori  and  a 
posteriori  error  estimates  of  finite  element  solutions  of  (1.2)  on  regular  branches  and  on  branches 
around  turning  points  in  the  case  that  the  number  of  parameters  is  one,  that  is,  A  C  R.  Since 
our  formulation  of  parametrized  nonlinear  equations  includes  strongly  nonUear  equations,  our 
theory  is  an  essential  extension  of  the  prior  works. 

In  this  paper  we  present  the  theory  of  a  priori  error  estimates.  In  [TBl]  the  theory  of  a 
posteriori  error  estimates  and  several  numerical  examples  will  be  given.  In  the  following  the 
outline  of  this  paper  is  described. 

First,  we  show  that  the  exact  and  finite  element  solutions  of  (1.2)  form  one-dimensional 
smooth  manifolds.  If  F  is  mildly  nonlinear,  showing  that  solutions  form  manifolds  would  not 
be  very  difficult.  If  F  is  strongly  nonlinear,  however,  it  would  become  very  difficult,  or  F  would 
not  be  even  differentiable  in  A  x  Hq{J). 

Therefore,  we  redefine  (1.2)  and  (1.3)  using  the  Sobolev  space  Then,  F  becomes 

as  smooth  as  the  functions  a  and  /,  and  it  is  a  Fredholm  operator  in  a  certain  open  set.  From  the 
Fink-Rheinboldt  theory,  we  conclude  that  the  exact  and  finite  element  solutions  form  smooth 
manifolds  under  suitable  conditions. 

Next,  we  prove  several  a  priori  estimates  of  finite  element  solution  manifolds  of  (1.2)  using  the 
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extended  implicit  function  theorem  due  to  Brezzi,  Rappaz,  and  Raviart  [BRRlj.  As  mentioned 
before,  we  need  to  take  the  Sobolev  space  as  the  stage  of  the  error  analysis  of  finite 

element  solution  manifolds.  However,  using  J)  in  the  formulation  make  the  finite  element 

analysb  difficult.  So  we  have  to  come  up  with  several  new  tricks  to  overcome  this  difficulty.  The 
following  is  the  most  essential  trick: 

Since  our  operator  F  is  defined  on  its  Frechet  derivative  is  a  linear  operator 

on  However,  D,^F  can  be  extended  to  an  element  of  C{Hl,  H~^)  and  thus  the  usual 

theory  of  finite  element  can  be  applied  to  D^F. 

Another  new  idea  is  ‘rotation’  or  ‘pivoting’  of  the  coordinate  to  handle  turning  points.  In 
[BRR2],  a  slightly  different  formulation  from  that  of  [BRRl]  was  used  to  deal  with  turning 
points.  In  Fink-Rheinboldt’s  theory,  certain  isomorphisms  were  introduced  in  the  formulation 
so  that  both  regular  branches  and  branches  around  turning  points  were  treated  simultaneously. 
In  this  paper,  we  put  an  auxiliary  equation  in  the  original  equation  (1.2)  so  that  the  enlarged 
operator  is  an  isomorphism  between  Banach  spaces  around  turning  points  or  on  'steep  slope’. 
Then  we  do  the  same  thing  what  we  do  on  regular  branches  to  the  extended  operator. 

In  this  paper  one-dimensional  case  is  dbcussed.  Under  certain  assumptions  the  results 
obtained  here  will  be  extended  to  two-dimensional  case  in  [TB2]. 

This  paper  is  a  revision  of  a  part  of  one  of  the  authors’  Ph.D.  dbsertation  [T]. 

2.  Preliminary. 

In  this  section  we  prepare  notation  and  a  necessary  lemma. 

Let  J  :=  (6,c)  C  R  be  a  bounded  interval.  For  a  positive  integer  m  and  a  real  p  €  [l,oo]. 
we  denote  by  the  usual  L^-Sobolev  space  of  order  m,  that  is, 

:=  {u  €  U'iJ) \d^u  €  L^{J),  0  <  k  <  m}  . 

We  define  the  norm  of  by 

k=0 

For  p  €  [l,oo],  we  define  the  closed  subspace  by 

:=  {«  €  |u  =  0  on  ay  }  . 

As  usual,  we  denote  J)  and  Wo’^(y)  by  and  Hq(J),  respectively. 
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Note  that  C^{J),  the  set  of  infiaitely  many  times  diifeientiable  (unctions  with  compact 

1  *1  -TO 

supports,  is  dense  in  Wq{J)  for  p,  1  <  p  <  oo,  but  if  p  =  00,  C^(J)  is  not  dense  in  Wq’  (J). 
By  the  Poincare  inequality,  the  norm 

(2.1)  •“  il“  lU' 

is  equivalent  to  the  norm  ||  •  ||w>.»  We  always  take  the  norm  (2.1)  for  in 

this  paper. 

For  1  <  g  <  oo  and  p  with  ^  =  1,  let  be  the  dual  space  of  with  the 

norm 

\\F\\w.u,  ;=  sup  I  p<  F,  I  >,  I,  Fe 

where  p<  •,•  >,  is  the  duality  paring  between  and  iyo’'’(J).  Then  we  have 

Lemma  2.1.  For  any  F  €  W~^'^{J)  with  1  <  p  <  oo,  there  exists  a  unique  u  e 
so  thdt 

p<  F,u>,=  y^uVdi,  Vv  6  Wo^’''(7).  □ 

Lemma  2.1  is  a  direct  consequence  of  [B,  Proposition  VIIL13]. 

In  notation  of  this  paper,  we  omit  ‘(/)’  from  the  notation  of  Sobolev  spaces  when  there  is 
no  danger  of  confusion.  Also,  we  write  <  •,  •  >  instead  of  p<  •,•  >q  when  the  setting  of  the 
duality  paring  is  obvious. 

Subscripts  like  a,  and  stand  for  partial  derivatives  with  respect  to  r  and  A,  respectively. 

3.  Formulation  of  the  Problem. 

In  this  section  we  formulate  our  problem  rigorously.  To  do  this  we  define  the  nonlinear  operator 
F  :  A  X  -  VF-i'*  by,  for  A  €  A  C  R  and  u  €  Wj'*, 

(3.1)  <  F(A,u),v  >:=  Jja(A,x,u'(x))v\z)  +  /(A,x,u(x))v(z)]c{x,  Vv  € 

where  <  •,  •  >  is  the  duality  pairing  between  and 

For  F  being  well-defined  and  smooth  we  require  several  conditions  for  a  and  /. 

A  function  t(':Ax7xR-*Ris  called  Carath^odory  continuous  if  rp  satisfies  the 
following  conditions:  for  (A,x,y)  €  A  x  7  x  R, 
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{ip{\,i,y)  is  continuous  with  respect  to  A  and  y  for  almost  all  x, 
rp{\,x,y)  is  Lebesgue  measurable  with  respect  to  z  for  all  A  and  y. 

If  ^(A,i,y)  is  Caratheodory  continuous,  ^(A,z,u(r))  is  Lebesgue  measurable  with  respect 
to  X  for  any  Lebesgue  measurable  function  u. 

Let  a  =  (01,02)  be  usual  multiple  index  with  respect  to  A  and  y.  That  is,  for  o  =  (01,02), 
D°a(X,x,y)  means  a(A,z,y). 

Let  d  >  1  be  an  integer.  For  o,  jo]  <  d,  we  define  the  maps  A“(A,u)  and  F°'(A,u)  for 
(A,ti)€  by 

(3.2)  A‘*(A,«)(z)  :=  D‘*a(A,z,u'(i)), 

(3.3)  F“(A,u)(i)  ;=  D“/(A,x,u(i)). 

We  then  assume  that 

Assumption  3.1.  For  all  a,  |o|  <  d,  we  suppose  that 

(1)  For  almost  all  x  ^  J,  D°‘a{\,x,y)  and  D^f{X,x,y)  exist  at  all  (A,  y)  6  A  x  R,  and  they  are 
Caiatbeodory  continuous. 

(2)  The  mapping  A**  defined  by  (3.2)  is  a  continuous  operator  from  A  x  IFo’*  to  F*,  and  the 
image  C  F®’  of  any  bounded  subset  U  C  A  x  Wq’°°  is  bounded. 

(3)  The  mapping  F®  defined  by  (3.3)  is  a  continuous  operator  from  A  x  Wq'°°  to  .  and  the 
image  F®({7)  C  of  any  bounded  subset  U  C  Ax  Wj’'*  is  bounded.  □ 

Assumption  3.1  is  satisfied  ifa,  /  :Ax7xR— ►R  are,  for  instance,  functions.  By  simple 

computation  we  obtain  the  following  Lemma. 

Lemma  3.2.  Suppose  that  a  and  f  satisfy  Assumption  3.1.  Then  F  defined  by  (3.1)  is  a 
C*  mapping.  Its  Frechet  derivatives  are  written  as 

<  Z?,F(A,«)^,t/ >  =  Jja,(X,x,u'(x))i/'v' +  fy(X,x,u(x))ri/v]dx, 

<  DxF(X,u)i],v  >  =  V  Jla\iX,x,u'ix))v' +  fx{X,x,u(x))v]dx, 

for  if  €  ^*td  r;  6  R.  Moreover,  we  have  the  following  estimates 

||Z?.F(A,«)(|^j^yj,oo  <  ||a,(A,z,«'(z))||i«.  +  ||/»(A,  z,  u(z))||i, , 

||Z?>F(A,u)||w'-i.oo  <  ||ai(A,z,u'(z))lk*  +  ||/A('^,2^,«(2;))||r,i-  C] 
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Now,  we  define  our  problem. 


Problem  3.3.  Undet  Assumption  3.1  with  d  >  1,  solve  the  following  equation:  Find  A  6  A 
and  u  6  such  that 

<  F(A,u),v>=0,  Vv€lVo^'\ 

where  F  is  defined  by  (3.1).  □ 


4.  Fredholm  Operator  and  the  Solution  Manifold. 


In  this  section  we  prove  that,  if  a(A,i,y)  satisfies  certain  conditions,  F  will  be  a  nonlinear 
Fredholm  operator  and  solutions  of  Problem  3.3  form  a  one-dimensional  differential  manifold. 
The  following  lemmas  are  essential. 

Let  p  €  (l,oo]  and  a  €  L°°.  Define  A  :  IVq’’’  — ►  W~^’P  by 

(4.1)  <  Au,v  >:=  J  a{x)u'{x)v'{x)dx,  Vu  € 

where  ^  ^  =  1,  and  <  •,  •  >  is  the  duality  pairing  between  and  Wq’’’.  Then  we  have 


Lemma  4.1.  Suppose  a~^  €  and 
and 


L 


dx 

o(x) 


96  0. 


Then  A  is  an  isomorphism  between 


f  dx 

Proof.  First,  we  prove  that  A  is  onto.  Let  cq  :=  /  —tt-  Take  an  arbitrarv  F  £  p.  Bv 

Jj  or(x) 

Lemma  2.1  we  know  that  there  exists  a  unique  ijj  €  Wq’^  such  that 

<  F,v  >=  j  ilf'{x)v'{x)dx,  Vu  €  Wq'’’. 

Let  Cl  := - f  ~y  ^ dx  and  u(x)  :=  /  dt.  Then  it  follows  that  u  6  VF}'’’,  and 

CO  Jj  a{x)  h  a(t)  ° 

<  i4u,  V  >=<  JP,  V  >  for  all  w  e  fVo'*'.  Hence,  A  is  onto. 

Next,  we  show  that  A  is  one-to-one.  Suppose  that  ui,ii2  €  and 

J  a{x)u'i{x)v'{x)dx  =  j  a{x)v,2{x)v'(x)dx,  Vv  £  VFq’^'. 

By  (B,LemmaVIII.l.],  there  is  a  constant  C2  such  that  o(i)(ui(i)  -  U2(x))  =  C2  for  almost  all 
X  £  J.  Since 

0  =  f  (u'(i)  -  «^(x))dr  =  C2  /  =  C0C2, 

Jj  Jj  a{x) 
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and  Co  ^  0,  we  conclude  that  u\{x)  =  ti2(^)  ~  “2-  That  is,  A  is  one-to-one. 

Since  A  is  continuous  and  bijective,  A~^  is  also  bounded  by  the  closed  graph  theorem. 
Therefore,  A  is  an  isomorphism  between  IVq'’*  and  □ 

Lemma  4.2.  Suppose  that  a~^  €  1“  and  j  =  0.  Then 

(1)  dimKerj4  =  1  and  KeiA  —  {<p  €  IVq’’^1<p'(x)  =  coa(r)~\  co  €  R}, 

(2)  Imyl  C  W~^'^  is  closed, 

(3)  dimCokeri4  =  1. 


Proof.  (1)  Let  <p(i)  :=  /  -7-r.  Then,  by  the  assumption,  we  have  C  and 

Jh  ot{t) 

<  A(fi,v  >=  0  for  all  u  €  Therefore,  tp  €  Ker.4. 

Conversely,  for  any  u  €  KerA,  there  is  a  constant  cq  such  that  a(x)u'(i)  =  cq  for  almost  all 
X  €  J.  This  implies  that  u  =  CQ<fi.  Hence,  (1)  is  proved. 

(2)  First,  we  define  the  subset  X  C  by 


Clearly,  X  is  a  closed  subspace  of  Wq'^. 

Let  T  €  £(1Vq’^,  be  the  bomorphism  defined  by  <Tu,v  >:=  J  u'v'dx,  'iv  € 

-  /■*  ‘4f'{t)  1  „ 

Let  X  :=  T{X).  Take  any  1/'  €  A",  and  define  u(i)  :=  /  .  'dt.  Then,  we  have  u  G  ,  and 

Jb  Offt) 

<  i4u,  V  >=<  Tij),  V  >,  for  all  u  €  Hence,  we  have  that  Imi4  D  X. 

Now,  take  any  rj  €  Wq’’'  and  define  7  by  7(1)  :=  j  {a{t)Ti'{t)  -  ci)dt,  where  cj  := 
aiii'dxl\J\.  We  check  that  7  €  Wq^  and  T7  =  A'q.  Moreover,  we  have  7  €  A  because 


/ 


Jj  a(i)  Jj  Jj  0(1) 


Hence,  we  conclude  that  Im.<4  =  A  and  Imi4  b  closed. 

(3)  As  before,  define  €  Wq'^  by  V'o(ar)  :=  j  a{t)~^dt.  Since  j {■tpQ/a)dx  =  j  a~^dx  ^ 
0,  we  have  V'o  ^  A.  Let  C2  :=  j  a~^dx  >  0.  Take  any  ^  6  Wj’**  and  let  C3  :=  J {rl''/Q)dx. 
Then  t/'  -  (c3/c2)t/'o  €  A  because 


L 


Ip'ix)  -  {C3/c2)tlfo{x) 


J,  =  /  -  s  ^ 

Jj  0(1)  C2  Jj 


dx 


=  0. 


Of(i)  Jj  0(1)  C2  Jj  a(i)2 

This  implies  that  for  any  ip  €  there  exbt  C4  €  R  and  ipi  ^  X  such  that  ip  =  ctipo  +  xL'\. 
The  uniqueness  of  such  decomposition  is  obtained  by  a  simple  computation. 
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Therefore,  we  showed  that  =  Imj4  ©span{T^o},  and  (3)  is  proved.  □ 

From  Lemma  4.1  and  Lemma  4.2  and  the  definition  of  Fredholm  operators,  we  finally  obtain 


Theorem  4.3.  If  a~^  €  L°°,  then  the  Uneu  opentor  A  defined  by  (4.1)  is  a  Fredholm 
opeiitor  And  ind^l,  the  index  of  A,  is  0.  □ 

Let  ns  now  return  to  our  main  problem.  We  define  the  subset  5  C  A  x  Wq'°°  by 
(4.2)  5:={(A,ti)€Ax<’'*la,(A,i,«'(i))-i6i;“}. 

Since  the  mapping  A  x  Wq'°°  3  (A,u)  aj,(A,i,ti'(i))  €  1“  is  continuous,  we  have 

Lemma  4.4.  If  a  And  f  satisfy  Assumption  3.1  with  d  >  1,  S  is  An  open  set  in  A  x 


Now,  from  the  standard  theory  of  Fredholm  operators,  we  obtain  the  following  theorem: 

Theorem  4.5.  Suppose  that  a  And  f  satisfy  Assumption  3.1  with  d  >  1.  Then  in  S,  the 
opeTAtor  F  :  S  —»  defined  by  (3.1)  is  a  nonlineAr  Fredholm  operAtor  of  index  1. 

Proof.  From  Lemma  3.2  and  Theorem  4.3,  the  operator  DuF{X,u)  :  Wq’°°  —  jy-i  *  jg  Fred¬ 
holm  and  its  index  is  0  for  (A,u)  €  S.  Since  DF{X,u)  :  R  x  Wq’°°  —  is  written  as 

DF{X,u)(r},‘tp)  —  D^F{X,u)iIj  +  r)DxF{X,u)  for  n  €  R  and  ip  £  Theorem  4.5  is  con¬ 

cluded.  □ 

We  define  the  subset  1l(F,S)  C  S  by 

(4.3)  Jl(F,  S)  :=  {(A, «)  €  S\DFiX,  u)  is  onto), 
and  have 

Lemma  4.6.  For  Any  (A,«)  €  11{F,S),  dim  KerDuF’(A,  u)  is  at  most  1. 

Proof.  Assume  that  N  :=  dim  KerI?«F(A,«)  >  2  for  some  (A,u)  6  1l(F,S).  Note  that  the 
elements  of  KerI?F(A,  u)  Ate  solutions  of  the  linearized  equation 

(4.4)  DF{X,u){p,iP)  =  pDxF(X,u)  + D,FiX,u)iP  =  0,  p  £  R,  ip  £  . 
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If  DxF{X,u)  ^  ImZ)uf’(A,  u),  by  (4.4),  we  obtain 


KeiDF{X,u)  =  {(0,0)  6  A  X  0  €  KeiD^F{X,u)} 

and  dim  KerZ?f(A,  u)  =  N  >  1.  This  contiadicts  to  (4.3)  and  indf=l. 

Therefore,  we  should  conclude  that  DxF(X,u)  6  ImZ?tt/’(A,u).  Let  tpo  G  14^0’°°  be  such  that 
Di(A,u)  =  D^F{X,u)ipQ.  Then  we  obtain 

(4.5)  KerZ?/’(A,ti)  =  {(/^,-m0'o  +  0)€Ax  iyo*-“|/z€R,0€  KerZ>„/’(A,u)}, 

and  hence  dim  KerJ5F(A,  u)  =  AT  +  i  >  1.  Therefore,  we  get  a  contradiction  again,  and 
Lemma  4.6  is  proved.  □ 

The  elements  of  'Jl(F,S)  are  caJled  regular  points.  The  elements  of  F{H(F,S))  are  called 
regular  values. 

By  Theorem  4.5,  we  can  apply  the  Fink-Rheinboldt  theory  ([FRl],[FR2],[R])  to  the  operator 
F  and  obtain  the  main  theorem  of  this  section. 

Theorem  4.7.  Suppose  that  a  and  f  satisfy  Assumption  3.1  with  d  >  1.  Let 
e  6  F('fl(F,S)).  Then 

M  =  M^  :=  {(A,  u)  €  7l(T,<S)  |  F(X,  u)  =  e} 

is  a  one-dimensional  C^-manifold  without  boundary.  Moreover,  for  each  ( A,  u)  G  M,  the  tangent 
space  T(x,u)M  at  (A,u)  is  KetDF(X,u). 

Therefore,  if  0  G  F{Tl{F,S)),  the  solutions  of  Problem  3.3  form  a  one-dimensional  C^- 
manifold  without  boundary  in  3Z(F,S).  □ 

In  the  sequel  of  this  paper  we  always  assume  that  0  G  F(7Z(F,S)). 

Now,  let  us  consider  the  linearized  equation  (4.4).  From  Lemma  4.6,  we  would  have  four 
cases  for  (A,u)  G  7l(F,S). 

Case  1.  KetDuF{X,u)  =  {0}  and  DxF{X,u)  G  ImZ?„F(A,«). 

In  this  case,  by  the  implicit  function  theorem,  there  exists  a  unique  map  A  9  A  u(A)  G 
Wq'  such  that  F(A,«(A))  =  0  for  any  A.  Hence,  this  case  corresponds  to  regular  branches. 
□ 

Case  2.  dim  KerZ?uF(A,  u)  =  1  and  DxF{X,u)  ^  lmDuF{X,u). 
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In  Case  2,  using  the  well-known  Liapunov-Schmidt  reduction  (see,  for  instance,  [GS]), 
we  can  show  that  this  case  corresponds  to  (general)  turning  points.  □ 

Case  3.  Ker£)u^('^. «)  =  {0}  and  DxF{X,u)  ^  ImZ)uf’(A, u). 

By  a  similar  argument  to  the  proof  of  Lemma  4.6,  we  see  that  this  case  cannot  happen.  □ 
Case  4.  dim  KerZ)a/’(A,  u)  =  1  and  DxF{X,u)  €  lmDuF{X,u). 

By  (4.5),  we  have  dimKerDJ’(A,  ti)  =  2  and  dimCoker2>F(A,  u)  =  1.  Hence,  this  is  not  the 
case  for  (A,u)  €  1l(F,S).  In  this  case  we  may  have  a  bifurcation  phenomenon.  □ 

By  the  above  consideration  we  now  know  that 

{X,u)^TliF,S)  44  (A,u)  €5  and  £)P(A,u)  €  AR  X  is  onto, 

=>  we  have  either  Case  1  or  Case  2. 

5.  Regularity  of  Solutions. 

In  this  section  we  examine  the  regularity  of  the  solutions  (A,u)  6  Mq.  To  do  it  we  need 
additional  assumptions.  Let  p',  2  <  p*  <  oo  be  taken  and  fixed. 

Assumption  5.1.  Under  Assumption  3.1  with  d>l,  we  assume  that 

(1)  For  all  X  €  A,  the  functions  a(A,-,  •),  aj,(A,-,  •) :  7  x  R  — ►  R  are  continuous. 

(2)  For  all  (A,y)  €  A  x  R,  there  exist  ax(X,x,y)  for  almost  all  x  £  J  and  are  Caratheodory 
continuous. 

(3)  The  composition  functions  /(A,  x,  u(x)),  ax(X,  x,  u'{x))  are  in  L^'  for  any  (A,  u)  e  A  x  . 
Moreover,  for  any  bounded  subsets  K  C  Ax 

{/(A,i,«(i))  €  (A,u)  €  A'},  {o,(A,z,«'(r))  €  L^'\  (A,u)  6  A'} 

are  bounded  in  U’ .  □ 

Lemma  5.2.  Let  (A,u)  6  Mq.  Suppose  that  Assumption  3.1  and  5.1  hold.  Then  u  € 

C\J). 

Proof.  Define  /o  by  /o(x)  :=  — /(A,  i,  u(i)).  By  Assumption  5.1(3),  we  have  fo  6  L^' .  Now. 
consider  the  following  equation 

(5.1)  ^'{x)v'(x)dx  =  /o(r)u(i)di,  Vu  €  /fo¬ 


il 


There  exists  a  unique  solution  $  6  of  (5.1).  Thus,  we  have 

(5.2)  a(\x,u\x))  = 

where  H{x)  :=  #'(i)  +  ci  with  some  constant  ci. 

Now,  for  a  fixed  A  6  A,  define  the  function  G  :  7  x  R  — ►  R  by  G{x,y)  :=  a{X,x,y)  -  H{x). 
Note  that  Gy{x,y)  =  ay(A,i,y)  and,  by  Assumption  5.1(1),  G  and  Gy  are  continuous.  Also 
we  remark  that,  for  almost  all  iq  €  7  and  yo  :=  tt'(io)*  we  have  G{xo,yo)  =  a{X,xo,u'{xo))  - 
H{xq)  =  0  and  ay(A,a:o,yo)  =  <iy{^,XQ,u'{xQ))  #  0  because  (A,u)  6  Ad  C  5  and  (4.2). 

Therefore,  by  the  implicit  function  theorem,  we  conclude  that,  in  the  each  neighborhood  of 
Xq,  there  exists  a  unique  continuous  function  T  such  that  T’(io)  =  “^(2^0)  ‘‘•nd 


G(x,T(x))  =  a(A,i,r(i))  -  Hix)  =  0. 

This  means  that  T(x)  =  u'(x).  Hence,  u'{x)  is  continuous  for  all  r  e  7.  □ 
The  following  is  the  main  theorem  of  this  section. 


Theorem  5.3.  Under  Assumption  3.1  and  5.1,  we  ha,ve  u  €  for  ail  (A.u)  e  Mq. 

Moreover,  for  all  bounded  closed  subsets  M  C  Mo,  there  exists  a  constant  K(M)  such  that 
sup_  ||u||vvrj,p.  <  Ji'(M). 

(\,u)€M 

Proof.  Let  H  €  be  defined  by  (5.2).  For  small  5  >  0  we  write 


ff(x  +  6)-ff(x) 
6 


ay(A,  x  +  6,  u\x  +  5)  +  €{u\x  +  5)  -  u'(r)))“ 

0 

a(X,x  +  6,u'(x))  -  a(X,x,u'(x)) 


with  0  <  e  <  1.  Since  If  €  ,  H'  exists  at  almost  all  z  €  7  and  H'  €  L^’  ■  From 

(A,  w)  €  Mo  C  S,  it  follows  that  a(A,  z  +  5,  «'(z  +  6))~^  €  L°°,  that  is, 


|a(A,z  +  6,u'{x  +  5))|  >  7  >  0  for  any  x  -k-  6  &  J. 

By  Lemma  5.2,  we  have  w'(z  +  5)  — ►  u'{x)  as  5  -*  0.  Hence  we  obtain 

(5.3)  |a»(A,  z  +  5,  u'(z  +  5)  +  e(u'(z  +  5)  -  u'(z)))|  >  0 
for  all  z  €  7  and  sufficiently  small  5  >  0,  and 

(5.4)  lim  ay(A,  z  +  6,  u'{x  +  S)  +  e(a'(z  +  5)  -  a'(z)))  =  ay(A,  z,  «'(z)) 
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because  of  Assumption  5.1(1).  By  (5.3),  (5.4)  and  Assumption  5.1(2),  we  conclude  that  ,  for 
almost  all  I  €  J,  liin(u'(i  +  5)  —  u'(i))/6  exists  and 

(5.5)  «"(i)  =  lim  +  <i)  -  =  tf'(x)-».(A.x,u'(r)) 

«-o  (5  ay(A,i,u'(i)) 

Since  H' ,ax{X,x,u'{x))  €  I*’*,  we  obtain  u"  €  i***  and  u  6 

Now,  let  A4  C  A4o  be  a  bounded  closed  subset.  Then,  we  have 

(5.6)  sup{||ay(A,i,u'(2:))“^||i,oa;(A,u)  e  M}  <  oo. 

The  last  part  of  Lemma  5.3  is  obtained  by  (5.5),  (5.6),  and  Assumption  5.1(3).  □ 

6.  Finite  Element  Solution  Manifold. 

Recall  that  we  ate  considering 

Problem  6.1.  Find  A  €  A  and  u  €  such  that 

(6.1)  <  F(A,ti),u  >=0,  Vv€W^’^.  □. 

0 

Naturally,  we  define  the  finite  element  solution  of  (6.1)  in  the  following  way.  Let  S/,  C 
C  be  a  finite  element  space.  The  space  of  piecewise  linear  functions  is  an  example 
of  S/,.  We  define  the  finite  element  solutions  of  Problem  6.1  by 

Problem  Find  A/,  €  A  and  u/,  €  such  that 

<  ■f’(‘^a>«a),Vfc  >=  0,  Vvfc  6  S/,.  □ 

Then,  using  the  Fink  and  Rheinboldt  theory,  we  will  show  that  the  solutions  of  Prob¬ 
lem  6.1 /<£  also  form  a  differential  manifold. 

Let  (•,•)  be  the  inner  product  of  defined  by  (u,u)  :=  /  u'v'dx  for  u,v  e  hfo-  Since 
SkC  H^,  we  define  the  canonical  projection  IIji,  :  hfj  —  S/,  by  (V'  -  Rhi>,Vh)  =  0,  Vu;,  6  5^  for 

We  see  the  following  equivalences.  Define  an  isomorphism  T  €  £(W~*  °°,  Wj'*)  by  < 
V,  V  >=  {Tt),  v),  Vu  €  Wq'^  for  ri  £  Then,  we  observe  that,  for  any  £  Sh  and  v  £ 

(6.2)  <  f{Xh,Uh),Vh  >=0  «=>  <  F(Afc,ri;,),  Ilat/ >=  0 

4=»  {TFiXh,Uk),llhv)  =  0 
4=>'  {EkTF{Xk,UH),v)  =  0 
<=»  <T~^llhTF{Xh,UH),v  >=0. 
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Since  Hq  b  dense  in  we  conclude  that  Problem  6. Iff;  is  equivalent  to 

O 

Problem  Find  A^,  €  A  and  Uf  €  such  that 

<  FH{Xh,tih),v>=Q,  Vu€^Vo’^ 

where  P*  :=  T-^UhT  6  and  F*(Afc,Ufc)  :=  ).  □ 

Our  formulation  of  Problem  seems  to  depend  on  T  and  Ilfc.  However,  we  claim  that, 
even  if  we  take  other  Ta,  and  define  the  finite  element  solutions  by  <  7’~^njTaP(Aft,  u/,),  v  > 
=  0  for  all  t/  €  this  formulation  is  equivalent  to  Problem 

Let  a  6  L°°  be  such  that  Qr(i)  >  e  >  0  for  all  x  6  7,  where  e  is  a  constant.  Let  (•,  •)„  be 
the  inner  product  of  Hq  defined  by  (u,  «)a  :=  j  au'v'dx  for  u,  w  €  Hq.  Define  the  isomorphism 

Ta  €  £(IV~^’‘^,PVq'°°)  by  <  Tj,v  >=  (Tai],v)a,  Vw  €  Wq’^  for  77  6  Also,  define  the 

0  • 

canonical  projection  :  Hq  by  (rp  —  =  0,  6  Sh  for  7p  6  Hq,  By  the  same 

manner  as  in  (6.2),  we  observe  that,  for  any  and  any  v  €  Hq, 

^  P(Afc,u/t),  Vh  >=  0  -i— r<  Tq  TaP(A^,  v^),  u  >=  0. 

Therefore,  with  the  definition  :=  T'^HgTa,  we  conclude  that 

(6.3)  PfcP(Afc,«fc)  =  0  <<=> /^P(Afc,Ufc)  =  0. 

Hence,  our  claim  is  demonstrated. 

We  will  see  that  these  observation  is  very  important  fc;  our  a  priori  error  estimates  because 
(6.3)  guarantees  that  we  can  take  any  a  €  I**  (that  is,  {Tq,  H^))  such  that  a  >  f  >  0  in  our 
error  analysis. 

In  the  statement  of  Problem  6.1J-£,  we  defined  Fh  :  A  x  —*  Following  the 

Fink-Rheinboldt  theory  we  extend  to  A  x  Wq'°°.  Define  pj  :  A  x  Wq'°°  —  by 

P?(A, «)  :=  (/  -  PS)Ta^u  +  P^F{X,  u), 

where  /  is  the  identity  of 

Lemma  6.2  ([R, Lemma  5.1]).  The  operator  7^  satisfies  the  [allowing: 

(1)  7X(A,u)  =  0  for  some  (A,u)  e  A  x  Hq  Hind  only  if{X,u)  €  A  x  5/,  a.nd  Fh(X,u)  =  0. 

(2)  TJ'  is  a  Ftedholm  opentor  of  index  1  on  A  x  Hq.  □ 
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By  Lemma  6.2,  we  have  the  following  theorem  as  a  consequence  of  the  Fink-Rheinboldt 
theory. 

Theorem  6.3.  Suppose  that  F  is  impping  (d  >  1).  Then  the  set  of  the  finite  elements 
solutions  of  Problem  6.  Ife> 

Mh  :=  {(Afc.ufc)  6  X  HI))\Fk{K,uk)  =  o}  , 

is  a.C^  minifold  without  boundary.  □ 

7.  A  Priori  Error  Estimates  of  the  FE  Solution  Manifold. 

Part  1:  Regular  Branches. 

We  ate  ready  to  start  to  consider  a  priori  error  estimates  of  the  FE  solution  manifold  Mh-  In 
the  consideration  of  error  estimates,  we  always  assume  the  following. 

Assumption  7.1.  We  assume  that 

(1)  Assumption  3.1  with  d  (i.e.  F  is  a  Fredholm  map). 

(2)  0  €  FiniF,S))  (i.e.  Xo  #  0). 

(3)  Assumption  5.1  (i.e.  u  €  ,  2  <  p*  <  oo  for  any  (A,  u)  €  Mo). 

(4)  Sk  is  regular  and  lim  inf  ||u  -  =  0,  for  any  u  e  ffo- 

o 

(5)  The  triangulation  of  S*  (in  one  dimensional  case,  the  partition  of  J  into  small  intervals) 
satisfies  the  inverse  assumption  [C,pl40].  □ 

In  the  sequel,  we  denote  by  Ilfc  :  Wq  — ►  S*  the  interpolant  projection.  We  also  denote  by 
C  or  C„  t  is  non- negative  integers,  generic  constants  which  are  independent  of  h  >0. 

The  main  tool  of  our  a  priori  error  estimates  is  the  following  implicit  function  theorem  due 
to  Btezzi,  Rappaz,  and  Raviart  [BRRl, Theorem  1]. 

Theorem  7.2.  Let  X,  Y  and  Z  be  Banach  spaces.  Let  S  C  X  and  y  :  S  —  F'  a  function 
defined  in  S.  Let  /  be  a  C*  mapping  defined  in  a  neighborhood  of  S  x  y{S).  Suppose  that  the 
function  5  3  z  y(z)  €  Y  satisfies  the  uniform  Lipschitz  condition;  there  exists  a  constant 
such  that 

||y(i)  -  y(z*)||y  <  Co||x  -  x‘\\x,  Vr,i*  €  5. 
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Suppose  in  addition  thit  the  following  hypotheses  hold: 

(i)  for  ill  xq  €  5,  Dyf{xQ,y(xQ))  is  in  isontorphism  ofY  onto  Z  with 


(ii)  we  hive 


sup  \\{Dyf{xQ,y{xo,)))~^\\c(z,Y)  ^ 
zo€5 


sup  \\D,f(x(i,y{xQ))\\c(x,Z)  <  C2, 

xo€5 

ind  there  exists  a  monotonially  incieising  function  Li  :  R'*'  — >  R'^  such  thit  for  ai/  zq  €  5  and 
ill  {x,y)  6  5(((zo,y(zo))) 


\\Df{x,y)  -  Df{xQ,y{xQ))\\c(XxY,Z)  <  ii(^)(ll*  -  lolU  +  llv  "  y(2;o)llK)- 

Then,  one  an  find  three  consUnts  a,b,d  >  0  depending  only  on  Cq,Ci,C2  ind  Li  such  thit, 
under  the  condition 


sup  ||/(zo,y(zo))||z  <  d, 

*o€5 

There  exists  a  unique  function  g  :  3a(xo)  — ►  Y  which  sitishes 

xo€S 


/(x,y(x))  =  0, 


ind  mips  i3(i(zo)  into  Bi,(y(xo))  for  xq  €  5.  Moreover,  we  hive  for  illxo  €  S  ind  ill  x  e  Ba(xo) 
||jf(z)  -  y(zo)||y  <  Ko{\\x  -  zolU  +  ll/(afo,  J/(io))||z). 
where  the  constant  A'o  >  0  depends  only  on  Ci.Cj.  □ 


Our  first  main  theorem  is  as  follows. 

Theorem  7.3.  Suppose  that  Assumption  7.1  holds  for  d  >  2.  Also  suppose  that,  at 
(Ao,  «o)  6  Afo»  D,tF{>>o,uo)  €  £(1^0’®’,  1^“^’°°)  is  in  isomorphism. 

Then,  for  sufficiently  smill  h  >  0,  there  ire  a  positive  >  0,  a  constant  6(Ao)  >  0,  ind  a 

-  • 

unique  map  [Aq  -  ca,,  Ao  +  ca,]  9  A  »-»  Ufc(A)  €  Sh  such  thit 

(7.1)  n(A,Ufc(A))  =  0, 

(7.2)  llufc(A)  -  nf.u(A)|l;,,  <  b(\o)hh\ 
for  any  Tf  with  0  <  q  <  j- 

Moreover,  we  have  for  some  consUnt  K{Xo)  >  0 

(7.3)  ||tifc(A)  -  «(A)||j5,.  <  /(r(Ao)llu(A)  -  n*«(A)||^. . 

Here,  the  constants  6(Ao)  and  A'(Ao)  are  independent  of  h  ind  A  €  [Aq  -  (Xo,Xq  + 
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Proof.  Since  the  proof  is  somewhat  long,  we  divide  it  into  several  steps. 

Step  1.  It  follows  from  Lemma  3.2  that,  for  (A,u)  €  Mq,  DF{X,u)  €  x 
and 

TAX  3  DF{\,u)  e  £(R  x  is 

I  Lipschitz  continaons  on  bounded  subset. 

We,  moreover,  claim  that,  if  D^F{\,  u)  6  ,  W"*’*)  is  an  isomorphism  for  (A,  u)  €  S, 

D^F{X,  u)  is  an  isomorphism  of  Hq  to  H~^  as  well. 

Define  Q,  R  e  CiH^,  H-^)  by 

2<  Q^,v  >2:— J^a{x)rij‘v'dx,  2<  A0,t;  >2:=  ^  £  Hq, 

where  a(i)  :=  a,(A,i,u'(i))  and  0{x)  :=  /y(A,i,«(i)). 

By  Theorem  4.3,  is  a  Fredholm  operator  of  index  0  and  R  is  compact.  Hence,  D^F{X,u)  = 
Q  +  iZ  6  is  a  Fredholm  operator  of  index  0.  Therefore,  if  KetDuF{X.u)  C  ^0 

trivial,  DnF{X,u)  €  H~^)  is  an  isomorphism. 

Let  Ip  £  Hq  he  such  that  DuF{X,u)ip  =  0.  This  means  that  -(af(x)^')'  +  3(x)tlj  —  0 
in  the  distributional  sense.  Since  P  £  XS  we  conclude  ip  £  by  a  standard  regular¬ 
ity  argument.  Hence,  ip  €  Wq’**  and  0  =  D^F{X,u)ip  £  Since  we  assumed  that 

DuF{X,u)  £  is  an  isomorphism,  we  obtain  ^  =  0.  Therefore,  our  claim  is 

proved,  o 

Step  2.  We  prepare  inequalities  which  we  will  use  later.  By  Assumption  7.1(5)  we  have 
the  inverse  inequality  [C,Theorem  3.2.6], 

(7.5)  l|vfc||^i,«  <  CaA-^llvfcil^,  <  e  5*, 

for  any  q  >  0. 

Since  D^F{Xq,  uq)  €  C{Wq'°°,  is  an  isomorphism,  we  conclude  by  the  implicit  func¬ 

tion  theorem  that  there  exist  Cj  >  0  and  a  unique  map 

(Aq  —  fj,  Ao  4-  €1)  9  A  t— ►  u(A)  6  Wg'°“ 

such  that  uo  =  u(Ao)  and  F’(A,u(A))  =  0.  Thus,  it  foUows  that 

(7.6)  ||nfc«(A-)  -  UHii{X)\\ai  <  C4IA-  -  A|, 

for  all  A,  A*  6  (A©  —  <1,  Aq  -b  €1 ).  In  (7.6)  we  used  the  fact  that  sup||n;k||rr»i  hI)  <  00  (see  e.g. 

h>0  ^  O’  0> 

[C, Theorem  3.1.6]). 
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By  Theorem  5.3  we  know  that  «(A)  6  and 

(7.7)  Cs  :=  sup  |l|«(A)||y^,j,^;  A  €  (Aq  -  y,  Aq  +  y]|  <  oo. 

Thus,  by  [C, Theorem  3.1.6],  we  see 

(7.8)  ||«(A)  -  nfcu(A)||^^..oo  < 

Let  ao(a?)  :=  Oy(AoiiJ.«o(i))-  Since  a,(Ao,x,«o(x))"^  6  ((Ao,«o)  €  Mo  C  S),  uo  € 

C  C’*[0, 1],  and  Assumption  5.1(1),  we  can  assume  without  loss  of  generality  that 

(7.9)  ay(Ao,  X,  Uq(x))  >  3^o  >  0- 

Then,  we  define  the  bilinear  form  A  :  Hq  x  R  by 

A(u,v)  :=  j  ao(x)u’v'dx,  Vu,i;  6  Hq. 

O 

Also,  we  define  the  canonical  projection  11°  ;  /Tg  by 

(7.10)  A(u-n2u,t;O  =  0, 

for  u  €  Then,  it  follows  from  Assumption  7.1(4)  and  (7.9)  that 

(7.11)  Um  ||«  -  njtiii*.  =  0, 

Now,  define  To  ;  by 

2<  V  >2=  ao(x){Toi)'v'dx,  Vu  6  Hq, 
for  ♦  €  and  define  7^  :  A  x  -*  by 

rJ(A,»)  :=  (/  - 

where  I  is  the  identity  map  of  and  :=  ro"^n2To  €  £(tV-^’“, 

Note  that,  by  Lemma  6.2,  7j(A,  u)  =  0  if  and  only  if  (A,u)  €  Mh- 
By  the  definition  of  11^  and  P^,  we  immediately  get 

(7.12)  C«  :=  sup||i^||£(jy-i  ,(,-1)  <  oo. 

a>o 

It  follows  from  (7.7),  (7.12),  and  [C,Theorem3.1.6]  that 

(7.13)  ||T;(A,nfc«(A))||^-.  <  Ce||/’(A,«(A))-F(A,nfc«(A))||ff-i 

~  (jo  “  nfcu(A)|(;f^ 

<  C7A,  VA  €  [Ao  -  y ,  Ao  +  y], 
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where  :=  D„/'(A,(1  -  t)ti(A)  +  tnfcu(A))  €  C(Hq,H~^).  Therefore,  we  obuin 

(7.14)  sup  {/i-<^+’')||7^(A,n*«(A))||^-.;A  6  (Ao  -  Ao  +  |-]J  <  I™  =  0, 

for  any  rj,  0  <  v  <  ^-  o 

Step  3.  We  claim  that  there  exist  a  positive  >  0  and  a  constant  Cg  >  0  independent  of 
h  >  0  and  A  €  [Ao  -  cj,  Aq  +  fa]  with 

(7.15)  ||D,T^(A,  nfc«(A))t;fc||//-i  >  ,  Vv^  6  S*. 

First,  we  note  that,  by  (7.4)  and  Step  1,  the  mapping 

(Ao  -  fi,  Ao  +  €i)  9  A  »-►  (Z?*/’(A, u(A)))“*  6  Hq) 

is  continuous.  Thus,  we  set 

Next,  we  write 

(7.16)  Z?.7;(A,ftAu(A))t;*  =  D^F{X,u{X)Wk' 

+i^(i>./’(A,fifcu(A))  -  i?./'(A,u(A)))vr 
-(/  -  /^)(-ro-'  +  D,F{Xo,uo))vh' 

+(/-  i^)( £)./’( Ao,uo)  -  D^F{X,u(X)))v["\ 

Let  us  examine  the  each  term  of  (7.16).  On  (a),  we  immediately  see 
(T  Ha)  ||X).i’(A,«(A))ufc||^-,  >  . 

On  (b),  it  follows  from  (7.4),  (7.8),  and  (7.12)  that 
(7.17b)  ||/^(i?./’(A,  nfcu(A))  -  Z?./’(A,  u{X)))vh\\„-i 

<  C,||n;,«(A)  -  u(A)||^..«||vfc||^.  <  CC5C9/i*-j!*-I|t;,||^.. 

On  (c),  we  remind  that 

2<  {-Tq^  +  DuF{Xo,uo))t{f,v  >2=  fy{Xo,x,uo{x))xi;vdx, 
and  —Tq^  +  D«F(Ao,uo)  6  C{Hq,  H~^)  is  compact.  Therefore,  we  conclude  that 
(7.17c)  lijn  11(7-  /^)(-ro-‘  +  7?.F(Ao.uo))||£(/r../,-.)  =  0 
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because  of  (7.11). 

On  (d)  we  observe  the  following.  By  (7.4)  there  exists  a  constant  C\q  such  that 

llO„f’(A'‘, u(A  ))  —  Oa/’(A, «(A)))|£(j^^  ^_i  j  <  CiolA*  —  A|. 

Take  ej  >  0  so  that  ej*  >  2wCiosup||/  -  P^Wau-^  H-')-  Then  we  have 

A>0 

(7.17d)  11(/  -  P?)(0.f(Ao.uo)  - 

for  any  A  6  [Aq  -  €2.  -^o  +  fj]- 

From  (7.16)  and  (7.17),  we  obtain 

llD.7t(A,nfcu(A))t/fcil^-.  >  ^:i:^-fi(h)^  IKII^,, 

with  lim  6(h)  =  0.  Therefore,  we  prove  the  claim  (7.15)  for  sufficiently  small  h.  o 
fc— 0 

Step  4  Again,  we  prepare  a  few  inequalities.  It  follows  from  (7.8)  that 

(7.18)  ||/)AT;(A,n*u(A))||^-.  <Cs(||i>iF’(A,u(A))||^..  +Cnh'-i^)  <Cj2, 

for  all  A  €  [Aq  ~  Aq  +  and  /i  >  0. 

By  (7.4)  and  (7.12)  we  have 

IP7^(A-.  v;)  -  Z?7;(A,n,u(A))||^,,^,_^..) 

<  CiadA*  -  A]  +  llt/J  -  ^^,u(A)l|^yJ,oo),  Vu;,  €  S^, 

where  C13  =  Ci3(|A|,  (A*|,  ||KjJ|(j^,^oo)  is  independent  of  h.  Hence,  by  (7.5),  there  exists  a 
monotonically  increasing  function  £1  :  R"*"  -►  R'*’  independent  of  h  such  that,  for  all  A,  A*  6 
(Ao  -  Ao  +  and  all  uJJ  €  with 

h-(Hv)(|A-  -  A|  +  lit;;  -  nA«(A)||^. )  <  e,  0  <  7,  <  i, 

we  have 

(7.19)  ||i?7;(A%t;;)-  Z?7;(A,n,u(A))||^,,^.,^-.) 

<  L,iOh-^^*'’\\X-  -  A|  +  lit;;  -  n;.u(A)||/,. ).  o 

Step  5  This  is  the  final  stage  of  the  proof.  By  (7.6),  (7.14),  (7.15),  (7.18)  and  (7.19)  we 
can  apply  Theorem  7.2  to  the  operator  in  the  following  situation; 
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A”  =  R  with  norm 

Y  =  Sh  C  Hq  with  norm  , 

Z  =  ShC  H~^  with  norm  , 

5  =  [Ao  -  €a,  ,  Ao  +  6a*]  with  6a*  :=  min  62) , 

y(A)  =  n;,u(A), 

where  7  :=  J  +  r/  for  any  7,  0  <  7  <  J.  Since  ||>lfc||£(XxK,z)  =  Phll£(Bx/r,',/r-‘)  aJl 

•  • 

A/t  €  £(R  X  5h,5h)  and  (7.14),  there  exist  a  constant  6(Ao)  >  0  independent  of  h  and  a  unique 
-  • 
function  [Ao  -  6a*,  Ao  +  6a*]  3  A  •-»  u/i(A)  6  S\  such  that 

(7.20)  7t(A,  iifc(A))  =  0  and  l|ni.ti(A)  -  u;,(A)||k  <  6(Ao), 
and  the  inequality 

(7.21)  ||tifc(A)  -  nkti(A)||;,.  <  Ci4|iT!(A,  nAtx(A))|l;,-. 

holds  for  all  A  €  [Ao  —  6a*,  Ao  +  6a*].  From  Lemma  6.2,  we  get  (7.1)  and  (7.2)  immediately  from 
(7.20).  Combining  (7.13)  and  (7.21),  we  obtain  (7.3)  and  complete  the  proof  of  Theorem  7.3.  □ 

Corollary  7.4.  Suppose  that  the  assumptions  of  Theorem  7.3  hold.  Then,  there  exists  a 
coastint  A’i(Ao)  >  0  independent  ofh>0  a.nd  A  6  [Ao  -  6a*  ,  Ao  +  6a*]  such  that 

(7.22)  ||ti(A)  -  tia(A)||j^,,»  <  Kt{Xo)h^  0  <  »?  <  1, 
for  any  A  €  [Aq  —  6a*,  Aq  +  6a*]  and  sufficiently  small  h  >  0. 

Proof.  By  (7.2)  and  the  inverse  inequality  [C,Theorem  3.2.6],  we  have 

l|iia(A)  -  na«(A)l|^y..«  <  b(Xo)h\ 

0 

for  all  A  €  [A©  —  6a*,  Ao  +  6a,].  It  follows  from  (7.8)  that 

||«(A)  -  «a(A)||^..«  <  l|ti(A)  -  nA«(A)|l^,.«  +  lltiA(A)  -  nAu(A)|l^..» 

<  CCih^~^  +  6(Ao)h’'  <  31^1  (Ao)*’.  □ 

Theorem  7.5.  Suppose  thit  Assumption  7.1  holds  for  d  >  2.  Also,  suppose  that  Mo  C 
Mq  is  i  compact  regular  branch,  that  is,  there  is  a  compact  interval  A  C  A  and  C-  map 
A  3  A  ►-*  u(A)  6  such  that 

Mo  ~  |(A,u(A))  €  Mq  I  /),/’( A,  u(A))  is  an  isomorphism  for  VA  €  a|  . 
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Then,  for  sufficiently  smaU  h  >  0,  there  exists  the  corresponding  finite  element  solution 
branch  Mh  C  Mh  which  is  parametrized  by  the  same  A  €  A  and 


(7.23) 

||nhu(A)-UH(A)t|^. 

< 

Koh2-^\ 

(7.24) 

||«(A)  -  «a(A)||;,^, 

< 

iirallu(A)  -  nA«(A)||^. 

(7.25) 

< 

K2h^ 

for  all  X  €  A,  tx(A)  6  Mo,  Ufc(A)  €  Mh,  and  j)  with  0  <  t)  <  ^.  Here,  A'o,A'i,A'2  >  0  are 
constants  independent  of  h  and  A. 

Moreover,  we  have 

(7.26)  MhCll{F,S). 

Proof.  From  Theorem  7.3  and  Corollary  7.4,  (7.23),  (7.24),  and  (7.25)  are  obtained  imme¬ 
diately. 

To  show  (7.26)  we  just  have  to  realize  that  Z)uF'(A,  u(A))  6  £(Wo’“,  is  an  isomor¬ 

phism  for  each  A  €  A  and 

D*f(A,UH(A))=:  D„f(A,u(A))+  Bh, 

where  Bh  :=  DuF{X,Uh{X))  -  DuF{X,u(X))  and  jy-i.co)  —  0  as  /t  —  0  because  of 

(7.25).  □ 

Remark  7.6.  We  can  rewrite  (7.23)  and  (7.25)  as 

||nfc«(A)  -  Uh{X)\\if,  <  AoA»-‘.  ||«(A)  -  u;.(A)||^.,»  <  A'zA*-', 

where  (  >  0  is  an  any  small  number.  In  linear  cases,  with  certain  assumptions  of  regularity  of 
solutions,  we  would  have  error  estimates  like 

PfcU  -  <  Ch^,  ||u  -  Uhll^roo  <  Ch. 

It  is  not  very  clear  whether  or  not  the  convergence  rates  of  (7.23)  and  (7.25)  are  optimal. 
We  might  be  able  to  improve  the  convergence  rates  with  further  assumptions  for  the  regularity 
of  solutions  (Assumption  5.1  might  not  be  enough  to  improve  (7.23)  and  (7.25)).  □ 

For  II  •  11^1  .^•-estimate,  we  have  the  following.  Suppose  that  we  have 

(7.27)  Um  llti  -  =0,  V«  6  Wj  "* , 
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where  11°  €  C(Wq^  ,  )  is  defined  by  (7.10).  For  example,  we  can  show  that  (7.27)  is  true 

for  piecewise  linear  elements. 

Theorem  7.7.  Suppose  that  assumptions  of  Theorem  7.5  a.nd  (7.27)  hold.  Then,  for  suf¬ 
ficiently  smaJl  h  >  0,  for  the  corresponding  finite  element  solution  branch  Mh  C  Mh,  we  have 
the  following  estimates: 


||n*u(A)-U4(A)||^.,,. 

< 

AohJ^^■"^ 

lKA)-MA)ll^^.,p. 

< 

KxMX)  -  n^u(A)|l^^i.,. 

IWA)  -  Uh(A)||^^..o<. 

< 

K2h\ 

for  ill  X  e  A,  u{X)  e  .V<o,  tifc(A)  €  Mh,  and  t]  with  0<7?<1-^.  □ 

8.  A  Priori  Error  Estimates  of  the  FE  Solution  Manifold. 

Part  2:  Ar.ound  Turning  Points. 

Let  ns  consider  a  priori  error  estimates  around  turning  points  and/oi  on  ‘steep  slopes’.  Basic 
idea  is  as  follows:  just  rotate  the  coordinate  ‘90-degrees’  and  do  the  exactly  same  thing  as  in 
Section  7. 

Recall  that  by  the  argument  in  Section  4  we  know  that  we  have  either  Case  1  or  Case  2  for 
(A,u)  6  .Vio  C 

Case  1:  KerD«F(A,u)=  {0}  and  DxF{X,u)  6  ImZ)BF(A,u). 

Case  2:  dimKeri?,F(A, «)  =  1  and  DxF{X,u)  ^  ImZ?ttF(A,u). 

Suppose  that  7  €  R  and  xo  €  J  are  given  in  a  certain  way  and  fixed.  Define  G  :  'R.[F,S)  — 
R  X  by 

G(A,u)  :=  (u(io)  -  7.  F(A,u)) 
for  (A,«)  €  7tiF,S).  Then,  we  have 

(81)  DG{X,u){p,rl>)  =  (i>{xo),fiDxF(X,u)  -I-  DuF(X,u)n;), 

for  M  €  R  and  if)  6  1^0’“.  First,  we  prepare  the  following  lemma. 

Lemma  8.1.  Let  (A,u)  €  '7J(F,5)  a,nd  (po>t('o)  the  basis  of  KeiDF(X,u).  Suppose  that 
DxF(X,u)  /  0.  Then,  we  have  ||ti'o||c«  >  0-  Moreover,  with  xq  €  J  such  that  il’oixo)  ^  0. 
DC{X,  u)  is  an  isomorphism  of  R  x  to  R  x  for  any  7  €  R. 
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Proof.  Suppose  that  we  are  in  Case  1.  It  follows  from  DxF(X,u)  ^  0  that 

••=  -D.F{X,u)-HDxF[X,u))  ^  0. 

Thus,  the  basis  of  Ket£>F(A,ti)  should  be  written  as  (Mo.Mot^'i).  Hence,  we  obtain  ll^ollc«  = 

ImoHIV'iHc®  >  0- 

Suppose  that  we  are  in  Case  2.  Let  ipQ  €  W'o  '*  be  such  that  KerZ?„f’(A, u)  =  span{^o}- 
Then,  (0,  V'o)  is  the  basis  of  KerZ)F(A,ti)  and  the  first  part  of  Lemma  8.1  is  trivial  in  this  case. 

Now,  let  us  consider  DG{X,u).  Let  (/i,  6  KerDG(A,ix),  and  xo  €  J  such  that  ip{xo}  #  0. 

Then,  there  exists  (5  €  R  such  that  (Mit/')  =  t^o)-  Thus,  it  follows  from  (8.1)  and  ipoixo)  ^  0 

that  =  (0,0).  Hence,  DG(A,«)  is  one-to-one. 

Let  (2,^)  be  any  element  of  R  x  Since  DF(A,u)  is  onto,  there  exists  € 

R  X  such  that  $  =  DF{X,  u){p,  0).  Hence,  we  obtain  (2,  $)  =  Z?G(  A,  u)((p,  0)  Vo)), 

where  6  :=  (z  -  0(io))/V'o(io)-  Therefore,  we  have  showed  that  DG{X,u)  is  onto  and  an 
isomorphism.  □ 

From  Lemma  8.1,  we  immediately  obtain  the  following  corollary. 

Corollary  8.2.  Let  (A,u)  6  11{F,S).  Suppose  that  DxF{X,u)  0.  Then,  lot  sufficiently 

0 

small  h  >  0,  there  exists  a  nodal  point  xq  €  J  oISh  suck  that  DG(X,u)  is  an  isomorphism.  □ 


Remark  8.3.  In  Lemma  8.2  we  showed  that  we  always  can  choose  a  nodal  point  of  Sh  so 
that  DG{X,u)  is  an  isomorphism  if  DxF{X,u)  ^  0.  For  example,  if  a  nodal  point  xq  ^  J  is 
taken  so  that  i>o{xo)  is  nearly  equal  to  |lt('ollc*>  DG{X,u)  is  an  isomorphism. 

Indeed,  the  manner  of  PITCON  of  choosing  the  continuation  index  is  consistent  to  the 
above  fact.  After  getting  a  point  (Ah,«j,)  €  Mh,  PITCON  computes  the  tangent  vector 
fk  =  (yo,"-,yfc)  €  ^1*®  solution  manifold  Mh  (remember  that  KeiD Fh(Xh,Uh)  = 

Tix,,  Then,  the  continuation  index  ic  is  taken  so  that  |y,^|  =  ||tA||oo-  In  our  case. 

O 

iyo,--,yk)  is  like  ipoh,i>0h{xi),  where  ii,...,ifc  are  nodal  points  of  5^,  (po/,,  vo/,)  is 

the  basis  of  KeiDFh{Xf^,utt),  and  y,^  =  0o#i(xo)-  (See  [R]  for  the  detail.)  Thus,  for  sufficiently 
small  h  >  0,  |V'oa(^o)|  would  be  very  close  to  ||t('o||co,  ^nd  0o(^o)  is  not  zero. 

Hence,  in  practical  computation,  we  may  expect  that  PITCON  takes  the  right  nodal  point 
zo,  »nd  DG{X,u)  £  £(R  x  x  IP"*’'*’)  is  an  isomorphism.  □ 
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Suppose  that  DxF{X,u)  ^  0  at  (A,  u)  6  Mo-  To  ‘rotate’  the  coordinate  we  define  the 
operator  H  :R  x  S  —  R  x  by 

Hi-y,  \,u)  :=  (u(io)  -  y,F{X,u)),  7  6  R,(A,u)  e  S, 

where  iq  €  J  is  taken  so  that  Z?(A^„)/f(7,  A,«)  =  DG(A,u)  e  C{R  x  x  is  an 

isomorphism. 

Note  that 

(8-2)  DHiy,  A,  u)(3,  t,  V')  =  {~s  +  ip(xo),  DF{\,  u)(t,  ^)) 

for  (3,  t)  €  R^  and  rp  €  ^o'°°-  Also,  note  that,  by  the  implicit  function  theorem,  for  each 
(A,u)  €  ^Vlo  such  that  DxF{X,u)  ^  0,  there  exist  cq  >  0  and  a  unique  map  {u(xo)  — 
60,  ti(xo)  +  fo)  9  7  (A(7), u(7))  €  Mq  such  that  (A,u)  =  (A(7o), u(7o))  with  70  :=  u{xq),  and 

■^(7)  A(7),  u(-i'))  =  (0,0),  that  is,  f(A(7),«(7))  =  0  and  u(7)(xo)  =  7o  for  any  7. 

Suppose  that  D\F{Xq,uq)  ^  0  at  (Ao,«o)  €  Mo-  Then  by  Corollary  8.2  there  e.xists  a  nodal 

O 

point  xo  €  J  of  5/,  such  that  f?(A,u)-^(7» Ao,tio)  =  DG(Xo,uo)  is  an  isomorphism  of  R  x  Wq'^ 
to  R  X  ly  !•<»  for  sufficiently  small  h  >  0. 

Theorem  8.4.  Suppose  that  Assumption  7.1  holds  lot  d  >  2.  Let  DxF(Xo,uo)  ^  0  at 

(Ao,«o)  €  Mo.  We  assume  without  loss  of  generality  that  there  exists  xq  €  7  such  that  xo  is  a 

0 

nodel  point  of  for  aii  sufficiently  small  h  >  0  a,nd  i?(A,u)'^f(7>  Aqi  uq)  is  an  isomorphism. 

Then,  there  exist  a  positive  co  >  0,  constants  b{Xo,uo),  k{Xo,uo)  >  0,  and  a  unique  C~  map 
[«o(xo)  -  fo,  Mo(xo)  +  eo]  9  7  •-*  (A/,(7),  fifc(7))  £  Ax  such  that 

(8-3)  FH(Xh(y),  «fc(7))  =  0, 

(8.4)  |A^(7)  -  A(7)|  +  |liiA(7)  -  n/.«(7)||//.  <  6(Ao,  uo)hh\ 

l^/i(7)  -  A(7)|  +  l|«fc(7)  -  «(7)||/f^  <  AT-'o,  Wo)||w(7)  -  n/,u(7)||Hj  - 

for  any  q  with  0  <  77  <  The  constants  6(Ao,  uq)  and  K{Xo,uo)  ate  independent  of  h  and 
7  €  [«o(xo)  -  fo,  «o(xo)  +  eo]- 

Proof.  The  manner  of  the  proof  is  exactly  same  to  that  of  Theorem  7.3.  We  divide  the  proof 
into  several  steps. 

Step  1.  It  follows  from  Lemma  3.2  and  (8.2)  that  D/f(y,X,u)  £  C{R‘  x  H^.R  x 

and 

(8  6)  I  ^  A  ^  *  9  (7,A,u)  ^  DH{y,X,u)  €  C{R^  x  H^R  x  R-^)  is 

1  Lipschitz  continuous  on  bounded  subset. 
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We,  moreover,  claim  that,  if  DxF(X,u)  ^  0  and  D(x,^)H{y,  X,u)  €  £(Rx  Wq’°°,Rx 
is  an  isomorphism  at  (7,  X,u)  €  R  x  1l(F,S),  X,u)  is  an  isomorphism  of  R  x  /Tq  to 

R  X  as  well. 

Let  H  :=  D(x,u)H{lA,u)  €  £(R  x  x  tf-i).  Then 

=  i'i’ixo),  DF{X,u){fi,^)),  fi  €R,i>  €  Hq. 

We  have  to  consider  two  cases. 

Suppose  that  we  are  in  Case  1,  that  is,  DuF{X,  u)  €  £(  is  an  isomorphism.  In 

this  case,  from  the  argument  of  Step  1  of  Theorem  7.3,  we  know  that  DaF{X,u)  6  C(Hq,  H~^) 
is  an  isomorphism.  Thus,  we  can  prove  our  claim  by  the  exactly  same  manner  of  the  proof  of 
Lemma  8.1. 

Next,  suppose  that  we  are  in  Case  2,  that  is,  dimKerQ  =  1  and  R  ^  ImQ,  where  Q  := 
DnF(X,u)  €  C{Wq°°  ,  and  R  ;=  DxF{X,u).  To  avoid  confusion  we  denote  D^F{X,u)  g 

We  first  show  that  KerQ  =  KerQ.  Obviously,  we  have  KetQ  C  KerQ.  Let  ^  6  KerQ  C  H^. 
Then,  we  have 

2<  Q'>(f,v  >2=  J [a(x)rp'v'  +  0(x)'tpv]dx  =  0,  Vu  €  Hq, 

where  0(1)  :=  ay(A,  i,  u'(i))  and  0{x)  ;=  /y(A,i,u(x)).  Hence,  by  a  simple  computation,  we 
conclude  that  ^  n  Hq  C  Hence,  KerQ  =  KerQ.  By  Theorem  4.3,  we  know  that 

indQ  =  0.  Thus,  we  obtain  dimCokerQ  =  1. 

Next,  we  want  to  show  that  DxF{X,u)  ^  ImQ.  If  DxF{X,u)  e  ImQ,  there  exists  some 
€  .^0 

J^a{x)ifj[v'dx  =  2<  DxF{X,u),v  >2  -  j  P{x)tpivdx,  Vu  e  ^o- 

Again,  by  a  simple  computation,  we  conclude  that  rpi  is  in  the  domain  of  Q.  This  is  a  contra¬ 
diction  because  we  assumed  DxF{X,u)  ^  ImQ. 

Since  we  showed  that  dim  KerQ  =  1  and  DxF{X,u)  ^  ImQ,  we  can  prove  our  claim  in  the 
same  way  as  in  the  proof  of  Lemma  8.1.  o 

Step  2.  We  prepare  several  inequalities  which  we  use  later. 

By  the  implicit  function  theorem,  there  exist  ci  and  a  unique  map 

(«o(io)  -  ei,«o(io)  +  «i)  3  7  •—  (A(7),«(7))  6  a  X  Wq'"* 
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such  that  (Ao,  uo)  =  (A(7o).  «(7o))  with  70.:=  wo(2;o).  «(7)(io)  =  7.  and  /’(A(7),  u(7))  =  0. 

As  in  Step  2  of  the  proof  of  Theorem  7.3,  we  easily  obtain  the  following  inequalities. 

(8.7)  |A(7’)  -  A(7)|  +  WflHuir)  -  nA«(7)||tf.  <  C15I7*  -  7|,  V7*,  7  6  h, 

where  /i  :=  [uo(afo)  -  M^o)  +  ^]- 

Let  00(2;)  :=  aj,(Ao,x,Uo(i)).  Again,  we  can  assume  without  loss  of  generality  that 

ay(Ao,x,tic(^))  >  3(5o  >  0. 

O 

Then,  we  define  the  canonical  projection  :  Hq  —<■  5a,  the  isomorphism  To  €  C(H~^ ,  Hq), 
and  the  map  7^  ;  A  x  as  in  Step  2  of  the  proof  of  Theorem  7.3.  Of  course, 

we  have  (7.11)  and  (7.12). 

Now,  define  :  R  x  5  —  R  x  by 

(8.8)  7^(7,  A,u)  :=(u(zo)-7,7t(A,u)),  7  6R,(A,u)e5. 

As  in  Step  2  of  the  proof  of  Theorem  7.3,  we  observe 

||:^(7,  A(7),  nAti(7))||«x/f-.  =  \fihu{y)ixo)  -  7|  +  ||7t (A(7),  nA«(7))||/f-. 

<  CisA,  V7e/i. 

Thus,  we  conclude  that,  for  any  77  with  0  <  t?  <  j, 

(8.9)  Un»  sup{/i~^7+’')||^(7,A(7),nAu(7))||Bx/f->}  =  0- 

Step  3.  We  claim  that  there  exist  a  positive  €2  >  0  and  a  constant  Cn  independent  of 
/»  >  0  and  7  6  [00(2:0)  ~  ^2. 00(2:0)  +  ^2]  such  that 

(8.10)  ||■D(A,n)^(7)  A(7).  nAo(7))(/2,OA)||Bxtf-'  >  C’lrdAtl  +  ||oA||/f|} ),  V^i  6  R.Vuh  G  Sa. 

First,  we  remark  that,  by  (8.6)  and  Step  1,  the  mapping 

(00(2:0)-  fi, 00(2:0) +  €i)  3  7*^  (-0(A,u)'*f(7,A(7),o(7)))~‘  €  £(R  X  H~\R  X  H^) 
is  continuous.  Thus,  we  set 

w  :=  max|I(Z)(j,„)/f(7,A(7),o(7)))"dl£(Bxtf-,«x//j)- 
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Next,  we  write 


(8.11)  Z?(A,„)^(7,A(7),nAu(7))(/i,i^h)  =  (vh(xo),  DF{Xiy),u('Y)){fi,Vk)) 

+(0,  (i?ft(A(7),  ny.r/(7))  "  DFiXh),  «(7)))(i^,  Vk))- 

On  the  first  term  of  right-hand  side  of  (8.11),  we  have 
(8-12)  ||(vh(io),  DF(X(y),  u(y))(tt,  Ufc))||Bx/f->  >  <*'“^1^1  +  )• 

On  the  second  term  of  right-side  of  (8.11),  we  write 

(8.13)  (i?7t(A(7),nA«(7))  -  Z)f(A(7),«(7)))(M,t^0 

=  -n(I-li)DxF(Xiy),uiy)f‘' 

+nP^(DxF{X(y),  n;.u(7))  -  DxFiXij),  u(7)))'‘' 
-»-Pfc°(Z?./’(A(7),  n,ri(7))  -  D^F{X(y),  uh)))v^’ 
-{I-P°){-To^  +  D^FiXo,uo))v["' 

MI  -  P^)iD.F{Xo,uo)  -  D,F{Xiy),u{y)))v'j;\  , 

Let  us  check  each  term  of  (8.13).  On  (a),  we  have 

(8.14a)  ||/i(/  -  P°)DxF{X{y),  u(7))|U-.  <  e(h)|/i|, 

with  lim  e(/i)  =  0  because  ||(/-  /^)i?A^(A(7),  u(7))||//-i  — ►  0  as  h  — •  0  uniformly  with  respect 
fc— *0 

to  7  on  /i  =  [uo(2;o)  ~  '<^o(io)  +  ^]. 

On  (b),  we  easily  obtain 

(8.14b)  ||/i/^(I>A/’(A(7),  n/.«(7))  -  DxF{Xiy),  u(7)))||h-  < 

On  (c)  and  (d),  we  immediately  get  (see  Step  3  of  the  proof  of  Theorem  7.3) 

(8.14c)  K(Z?„/’(A(7),  n;.u(7)  -  /?„F(A(7),  u(7)))t/fc||H-.  <  . 

(8.14d)  UrnJK/  -  P^)i-To^  -b  £>»F(Ao,uo))|la/f'.y/-)  =  0. 

On  (e),  by  (8.6),  there  exists  a  constant  C20  such  that 

||Z?«/’(A(7),«(7))  -  i?»i^(A(7*),u(7*))||ai/>,/f-)  <  C20I7  -  7*1, 

for  any  7,7*  €  Ii-  Take  63  >  0  so  that  >  2wC2oSup|I/  -  Ph\\c(H-'  //-')•  Then,  we  have 

fc>o 

(8.14e)  11(7  -  7»^)(D„f(Ao,«o)  -  .Dui^(A(7),«(7)))||£(/f-i,//->)  < 
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for  any  7  €  [uo(2;o)  -  €21  “0(2:0)  +  f2l- 
From  (8.12)  and  (8.14),  we  obtain 

IIAa.«)^(7, A(7),nfc“(7))(/2,“A)||«x/f-  >  (^  -  '5(/‘))(I/2|  +  ll“fcll//■), 
with  lim  6(h)  =  0.  Therefore,  we  prove  the  claim  (8.10)  for  sufficiently  small  h  >  0.  o 

A-»0 

Step  4.  Again,  we  prepare  a  few  inequalities.  By  (8.8),  we  see 

(8.15)  ||i?^:^(7,A(7),nfc“(7))||«xtf-  =  ll(-l,0)||,x/f->  =  1- 
Also,  we  immediately  obtain 

(8.16)  ||D:^(7',A-,u-)  -  i?:^(7,A(7),n,u(7))||£(B»x/r-,0lxif-.) 

<  C.iiir  -  7|  +  |A-  -  A(7)|  +  ||u*  -  n4u(7)llvv,..»). 

where  Cji  =  C’2i(|7*‘|,  |A*|,  |A|,  ||u•||^^^oo). 

Thus,  by  the  inverse  inequality  (7.5)  and  (8.16),  there  exists  a  monotonically  increasing 

function  X3  :  R'*’  -*  R'*'  independent  of  h  such  that,  for  all  7,7“  €  [“0(2:0)  -  “0(2^0)  +  ^]. 

« 

A*  €  A  and  vj  €  5),  with 

-  7|  +  |A“  -  A(7)|  +  IK  -  n/.“(7)||//.)  <  ^,  0  <  r,  <  1, 

we  have 

(8.17)  A-,  „•)  -  DHlii,  A(ir), 

<i;i(«)*-'^*'''(l7‘-7l  +  |A--A(7)|  +  ||v;-n»»(7)||a.).  « 

Step  5.  By  (8.7),  (8.9),  (8.10),  (8.15)  and  (8.17),  we  can  apply  Theorem  7.2  to  the  operator 
7^  in  the  following  situation; 

AT  =  R  with  norm  /»~*|7|, 

y  =  RxSfcCRx7o  with  norm  A“A(|A|  +  ||vfc||iiri), 

Z  =  Rx5fcCRx  with  norm  h'^^dAI  +  ||T,r^“fclU-' )» 

5  =  [uo(io)  -  €o,«o(2;o)  +  eo]  with  €0  :=  min(^,e2). 
y(7)  =  (A(7),  nfc“(7))- 
where  x  :=  ^  +  ^  for  any  tj,  0  <  tj  <  ^. 
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Since  ||^fcl|£(.vxK,2)  =  aJl  A/,  €  £(R^  x  S^.R  x  S^)  and  (8.9),  there 

exist  a  constant  b{Xo,  wo)  >  0,  independent  of  h,  and  a  unique  function  [uo(xo)  —  to,  uo(io)  + 
fo]  9  7  ^  (Afc(7).«fc(7))  €  R  X  Sfc  such  that 

(8.18)  ^(7,Afc(7),tifc(7))  =  (0,0),  1|(A(7),  nfcu(7))  -  (Aa(7), iifc(7))||y  <6(Ao,Uo), 

and  the  inequality 

(8.19)  |A*(7)  -  A(7)|  +  ||u/,(7)  -  nAu(7)IU,i  <  C22||^(7,  A(7),  n/iu(7))|l*x/r-'  • 

It  is  clear  that  (8.18)  implies  (8.3)  and  (8.4).  To  get  (8.5),  we  observe  that 

(8.20)  11^(7,  A(7),nkii(7))l|.x/f->  =  l|i^(^(A(7),  «(7))  - /’( A(7),  nAu(7)))||//-. 

<  c’23||«(7)  -  nfc«(7)ll/f^- 

Therefore,  combining  (8.19)  and  (8.20),  we  obtain  (8.5)  and  complete  the  proof.  □ 

By  the  same  way  as  in  Section  7,  we  obtain  the  following  propositions. 

Corollary  8.5.  Suppose  that  the  assumptions  of  Theorem  8.4  bold.  Then,  there  exists  a 
constant  A'i(Ao,uo)  >  0  independent  of  h  >  0  and  y  €  [uoC^^o)  ~  fo>  wo(io)  +  ^o]  such  that 

|A(7)  -  Afc(7)|  +  ||«(7)  -  «fc(7)|lnrj.«  <  hri(Ao,  uolh", 

for  any  0  <  q  <  j  and  y.  □ 

Theorem  8.6.  Suppose  that  Assumption  7.J  holds  for  d  >  2.  Let  Mq  C  Mq  be  a  connected 
compact  subset  with  the  following  properties: 

(1)  DxF{X,u)  ^  0  for  any  (A,ti)  6  A4o- 

(2)  There  exist  xo  €  J  such  that  Aa.  ,)/f(7,  A,u)  defined  by  (8.2)  is  an  isomorphism  for  all 
(A,  u)  €  Ado. 

Then  Mo  is  parametrized  by  y  =  u(io).  We  assume  without  loss  of  generality  that  the 
above  xo  is  a  nodal  point  of  Sh  for  all  sufficiently  small  h  >  0. 

Then  there  exists  the  corresponding  finite  element  solution  branch  Mh  C  M^  which  is 
parametrized  by  the  same  y\  that  is,  Ufc(7)(io)  =  7  and  A(A»,(7),  0^(7))  =  0  for  any  y. 
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Moreover,  we  hdve 


/ 


|a(7)  -  a*(7)|  +  l|n»«(7)  -  lik(7)ll».  < 

|A(7)  -  Afc(7)|  +  l|u(7)  -  tifc(7)||^^  <  /iL'3l|«(7)  -  nfcu(7)IU^  , 

|A(7)  -  Afc(7)|  +  l|u(7)  -  tik(7)||jj,^oo  <  K4h’>, 

Mh  C  1l(F,S), 

for  dll  y  =  ti(io),(A(7).«(7))  €  Xo,(A/k(7),«fc(7))  €  Mh,  and  ij  with  0  <  jj  <  ^.  Here, 
A*2.  Fj,  A\  are  positive  constdats  independent  of  h  dnd  y.  □ 

For  the  Wq^  -norm  estimate,  we  have  the  following  theorem  as  in  Section  7. 

Theorem  8.7.  Suppose  that  the  assumptions  of  Theorem  8.6  and  (7.27)  hold.  Then,  for 
the  corresponding  finite  element  solution  brdnch  we  have  the  following  estimates. 

|A(7)  -  Afc(7)|  +  ||na«(7)  -  tifc(7)|lvy..,-  ^ 

|A(7)  -  Afc(7)|  +  |1«(7)  -  Ufc(7)|lnr^^  <  A's||u(7)  -  n/,u(7)|lw-J.»>* - 
|A(7)  -  A„(7)|  +  ||«(7)  -  Ufc(7)||v^.,oo  <  K7h\ 

for  dll  y  =  u(io),(A(7),«(7))  €  Mo,  (Afc(7),  14^(7))  €  Mh,  and  ri  with  0  <  17  <  1  -  Here, 
R'i,  Ko,  A'r  are  positive  constants  independent  of  h  and  7.  □ 
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The  Laboratory  for  Numerical  Analysis  is  an  integral  part  of  the  Institute  for  Physical 
Science  and  Technology  of  the  University  of  Maryland,  under  the  general  administration  of  the 
Director,  Institute  for  Physical  Science  and  Technology.  It  has  the  following  goals: 

•  To  conduct  research  in  the  mathematical  theory  and  computational  implementation  of 
numerical  analysis  and  related  topics,  with  emphasis  on  the  numerical  treatment  of 
linear  and  nonlinear  differential  equations  and  problems  in  linear  and  nonlinear  algebra. 

•  To  help  bridge  gaps  between  computational  directions  in  engineering,  physics,  etc.,  and 
those  in  the  mathematical  community. 

•  To  provide  a  limited  consulting  service  in  all  areas  of  numerical  mathematics  to  the 
University  as  a  whole,  and  also  to  government  agencies  and  industries  in  the  State  of 
Maryland  and  the  Washington  Metropolitan  area. 

•  To  assist  with  the  education  of  numerical  analysts,  especially  at  the  postdoctoral  level, 
in  conjunction  with  the  Interdisciplinary  Applied  Mathematics  Program  and  the 
programs  of  the  Mathematics  and  Computer  Science  Departments.  This  includes  active 
collaboration  with  government  agencies  such  as  the  National  Institute  of  Standards  and 
Technology. 

•  To  be  an  international  center  of  study  and  research  for  foreign  students  in  numerical 
mathematics  who  are  supported  by  foreign  governments  or  exchange  agencies 
(Fulbright,  etc.). 

Further  information  may  be  obtained  from  Professor  I.  Babu§ka,Chairman,  Laboratory  for 
Numerical  Analysis,  Institute  for  Physical  Science  and  Technology,  University  of  Maryland,  College 
Park,  Maryland  20742-2431. 


